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ABSTRACT: 


An  interactive  computer  graphics  antenna  gain  pattern  computation 
and  display  program  for  real-world  antenna  systems  is  presented.   The 
use  of  the  program  as  a  teaching  tool  at  the  Naval  Postgraduate  School 
is  discussed.   Methods  for  applying  the  program  for  the  synthesis  and 
design  of  complex  antenna  systems  are  indicated.   Research  applications 
include  techniques  for  rapid  inspection  of  gain  equations  of  newly 
developed  antennas.   A  ship  motion  model  is  developed  for  studying  the 
effects  of  electrical  geometry  variations  caused  by  ship  motion  in 
heavy  seas  on  maritime  antenna  systems  and  a  dynamic  presentation  of 
pattern  variations  is  made. 
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A.   INTRODUCTION 

Anyone  who  has  attempted  to  correlate  the  actual  performance  of  HF 
and  VHF  wire  antennas  operating  in  a  real-world  environment  to  the  "highly 
sanitary,"  theoretical  radiation  patterns  and  gain  which  proliferate  in 
text  books  and  handbooks  will  immediately  recognize  the  need  for  a  simple 
method  of  predicting  antenna  performance  in  the  presence  of  the  earth.   This 
report  describes  the  development  and  use  of  a  computer  tool,  which  enables 
anyone  with  a  working  knowledge  of  Fortran  and  a  free-space  radiation  pattern 
and  mutual  impedance  formulation  to  analyze  and  design  antenna  systems  of 
arbitrary  orientation  above  a  specified  lossy  plane  earth. 

It  serves  as  both  a  teaching  aid  and  a  design  tool.   In  instructional  use, 
it  provides  quick  interaction  via  a  graphics  display  of  antenna  parameters 
and  radiation  patterns  plots.   In  seconds,  a  student  can  observe  the  perfor- 
mance of  several  popular  HF/VHF  antennas  in  any  plane  earth  configuration 
he  chooses.   Equivalent  digital  computer/plotter  turn-around  time  is  in 
excess  of  1  hour  and  manual  calculation  time  on  the  order  of  days.   The 
program  firmly  convinces  the  student  that  the  antenna  system  is  composed  of 
the  antenna  plus  its  environment. 

For  the  communications  system  designer,  rapid  evaluation  of  antenna 
systems  enables  him  to  choose  an  antenna  type  and  orientation  which  will 
enhance  the  performance  of  the  total  system  rather  than  arbitrarily  guessing 
which  antenna  package  might  be  suitable.   When  new  antenna  types  are  developed, 
their  radiation  pattern  and  impedance  equations  can  readily  be  added  to  the 
basic  calculation  package  and  the  full  potential  of  the  antenna  may  be  pain- 
lessly determined,  not  just  for  the  usual,  mystical  free  space  environment, 
but  for  the  surroundings  in  which  the  radiator  will  be  used.   Using  an  HF 
Ionospheric  Propagation  Prediction  program  which  can  return  optimum  radiation 
angles  for  a  specified  path  and  time,  a  systems  designer  can  synthesize  an 
optimum  antenna  pattern  for  the  particular  situation  under  investigation. 
The  optimum  may  be  manually  entered  and  each  design  iteration  compared  with 
the  optimum.   The  pattern  Save  and  Recall  options  are  used  to  arrive  at  the 
best  type  antenna  and  orientation  available  to  him. 

Previous  work  on  antenna  patterns  in  the  presence  of  ground  is  widely 
scattered  in  the  literature.   Specific  antenna  types  are  referenced  in  the 
appendices.   The  initial  incentive  for  this  investigation  was  to  increase 
the  usefulness  of  an  antenna  radiation  pattern  subroutine  for  HF  antennas, 
currently  in  use  as  part  of  a  HF  Ionospheric  Propagation  Prediction  program 
written  by  ESSAl.   The  gain  and  input  resistance  equations  from  this  report 
with  some  modifications  and  corrections  were  used.   Equations  programmed  are 
included  in  Appendix  C. 

The  program  as  presently  configured  assumes  a  current  distribution  on 
the  antenna.   The  consequences  of  this  are  small  errors  in  terminal  impedance 
with  a  corresponding  discrepancy  in  gain  value.   Radiation  patterns  are 
affected  very  little  by  the  differences  between  assumed  and  actual  currents. 
To  calculate  exact  current  distributions  would  be  prohibitive  in  both  time 
and  programming  effort,  considering  the  limited  worth  of  the  more  exact 
gain  figures  which  would  result. 


1  Ref:   ITS  78  Report. 


B.   BRIEF  DESCRIPTION  OF  PROGRAM 

The  program  consists  of  two  basic  parts: 

1.   The  solution  of  antenna  pattern  equations  for  radiators  of  arbitrary 
orientation  above  a  flat  earth  of  specified  ground  constants  (conductivity 
and  permitivity) .   Gain  values  are  calculated  from  terminal  impedance 
expressions  containing  self  (free  space)  and  mutual  (coupling  between' the 
antenna  and  its  image)  effects. 

2.   The  graphics  display  portion  which  displays  program  input  parameters 
specifying  antenna  type,  size,  orientation,  ground  constants  and  special 
features  such  as  recall  and  storage  of  patterns.   This  part  of  the  program 
also  generates  power  intensity  plots  vs.  azimuth  and  elevation  angles  and 
displays  these  radiation  patterns  on  the  graphics  screen.   Gain  values  can 
be  displayed  in  conjunction  with  the  patterns. 

Subroutines  for  special  functions  which  are  usually  found  in  wire  antenna 
patterns  and  impedance  formulae  and  numerical  quadrature  calculations  are 
included  for  the  convenience  of  persons  wishing  to  apply  their  own  specific 
antenna  to  the  program.   The  user  who  wishes  to  do  this  must  provide  pattern 
equations  for  his  antenna  for  arbitrary  orientation.   The  effect  of  ground 
reflections  and  the  selection  of  observation  angles  is  provided  by  the 
program  itself. 

Special  features  available  to  the  user  are: 

1.  Plotting  patterns  on  a  log  scale  vs.  linear. 

2.  Storage  and  recall  of  patterns  for  comparison  purposes. 

3.  Ability  to  generate  a  desired  pattern  shape  (via  light  pen)  which 
is  stored  and  recalled  for  comparison. 

4.  Simple  ship-ocean  model  for  dynamic  simulation  of  shipboard  antenna 
systems . 

When  the  program  is  used  at  USNPGS,  the  XDS  9300  digital  computer  core 
limitations  restrict  the  calculation  and  viewing  of  one  pair  of  cuts  in  the 
3  dimension  geometry  (i.e.  one  azimuth  rotation  at  one  specified  elevation 
angle  and  one  zenith  to  horizon  elevation  cut.)  Typical  time  for  a  pattern 
calculation  is  15  seconds,  with  a  simple  dipole  requiring  9  seconds  and  a 
vertical  monopole  with  ground  screen  up  to  2  minutes.   Seven  common  antennas 
are  currently  programmed: 

1.  Arbitrary  Tilted  Dipole 

2.  Vertical  Monopole 

3.  Vertical  Monopole  with  Ground  Screen 

4.  Inverted  L 

5.  Sloping  Long  Wire 


6.  Rhombic 

7.  Vertical  Half-Rhombic 

More  complex  programs  for  arrays  such  as  Yagis,  Log  Periodic  Dipoles 
and  Monopoles,  and  curtains  will  require  fairly  long  calculation  times  due 
to  the  extensive  mutual  impedance  calculations. 

C.   DYNAMIC  ANALYSIS  OF  SHIPBOARD  ANTENNAS 

The  ease  of  obtaining  the  effect  of  the  earth  on  antenna  performance 
prompted  the  investigation  of  the  programs  potential  to  display  the  effect 
of  typical  ship  motion  of  shipboard  HF  antenna  radiation.   The  equations  for 
an  arbitrary  tilted  dipole,  sloping  long  wire  and  vertical  monopole  were 
already  in  the  form  to  allow  variable  tilt  angle.   By  programming  a  ship- 
ocean  model  that  reorientates  the  antenna  with  ship  motion  (roll  and  pitch) , 
it  is  possible  to  show   slow  but  dynamic  pattern  changes  with  sea  surface 
as  a  function  of  sea  state  and  ship  direction  for  a  chosen  type  of  vessel. 
This  simplified  model  rocks  the  antenna  in  two  planes  as  the  ship  responds 
in  roll  and  pitch  to  ocean  waves  but  still  assumes  a  plane  ocean  reflecting 
surface.   For  medium  and  heavy  seas,  the  results  indicate  an  appreciable 
re-lobing  effect  and  show  that  the  variation  in  signal  at  a  particular 
observation  angle  may  be  as  high  as  20  db . 

This  is  an  additional  factor  which  should  be  considered  when  assigning 
locations  for  antennas  in  new  ship  designs.   Present  efforts  at  evaluation 
of  these  antenna  locations  do  not  include  sea  state  perturbations. 

The  next  stage  of  investigation  of  ship  motion  effects  will  include  a 
variable  geometry  for  the  sea  surface  to  replace  the  plane  shape.   At  low 
HF  the  effect  might  be  approximated  by  a  partially  random  "fuzzy"  surface 
while  for  UHF  the  distances  are  large  in  term  of  radio  wavelengths  and  the 
model  could  be  more  nearly  that  of  a  rolling  surface  contour. 

The  final  results  of  these  extensions  will  be  of  benefit  to  the  antenna 
locator,  as  previously  explained,  as  well  as  to  communications  managers. 
Depending  upon  the  sensitivity  of  the  total  communications  link  to  antenna 
lobe  structure,  communications  procedures  for  a  given  frequency  may  be  improved 
by  insight  obtained  from  the  investigation  of  sea  state  effects. 


D.   RECOMMENDATIONS 

The  radiation  equations  for  most  of  the  antennas  do  not  include  the 
arbitrary  geometry  factors  needed  for  ship  motion  effects  study  and  should 
be  expanded  to  include  them.   Gain  calculations  depend  upon  input  resistance 
which  in  the  case  of  the  sloping  long  wire  and  others  do  not  include  mutual 
effects.   Where  possible  and  where  warranted,  these  effects  should  be  included 
by  deriving  coupling  terms  for  input  impedance  equations.   (This  will  not 
alter  the  shape  of  the  radiation  patterns  and  affects  only  the  magnitude  of 
the  fields  and  hence  the  gain). 

Array  antenna  equations  exist  in  the  literature  and  should  be  carefully 
verified  and  adapted  for  inclusion  in  the  plotting  program. 


When  this  program  is  used  for  matching  patterns  produced  by  HF  propagation 
prediction  programs,  a  convenient  data  interfacing  technique  (such  as  tape) 
should  be  developed  for  use  between  the  graphics  system  and  the  lager  general 
purpose  digital  machine  used  in  the  prediction  calculations. 


APPENDIX  A 

DETAILED  PROGRAM  DESCRIPTION 
Appendix  A 

Section  III  contains  a  description  of  program  operation.   The  program 
is  divided  into  processors,  program  subsections  that  perform  the  major 
computational  tasks.   Processor  operation  and  interaction  are  described. 


1.   PROGRAM  OPERATION 

The  program  displays  a  data  and  option  command  input  format  at  the 
graphics  terminal  (see  figure  A-l) .   The  program  operator  enters  applicable 
parameters  values  for  antenna  geometry,  environment,  and  option  commands 
using  text  editing  techniques.   A  blank  graphics  block  is  then  displayed  at 
the  CRT.   Trial  patterns  may  be  manually  entered  in  this  block  using  graphics 
editing  techniques.   Manually  entered  patterns  will  be  displayed  with  all 
subsequently  computed  patterns  allowing  the  operator  to  compare  computed 
patterns  with  trial  patterns  on  the  CRT.   Exercising  the  reinitialization 
option  will  erase  the  trial  pattern. 

The  program  computes  horizontal  and  vertical  gain  patterns  and  displays 
them  at  the  graphics  terminal.   The  horizontal  pattern  is  computed  with 
zenith  constant  at  the  inputed  value  for  0'  and  azimuth  varied  from  1  to  360 
degrees  by  one  degree  increments.   The  vertical  pattern  is  computed  for 
azimuth  constant  at  the  inputed  value  of  $'and  zenith  varied  from  1  to  90 
degrees  by  1  degree  increments.   Linear  and  log  displays  are  available.   If 
a  log  display  is  not  ordered  with  the  log  pattern  option  command,  linear 
patterns  will  be  displayed. 

Patterns  are  saved  by  exercising  the  save  pattern  option.   Pattern 
vector  data  is  stored  in  the  digital  machine  in  a  save  array  when  save  is 
ordered.   Exercising  the  recall  option  will  cause  patterns  saved  in  the  save 
array  to  be  displayed. 

Use  of  save  and  recall  options  allows  simultaneous  display  of  saved  and 
current  patterns  for  comparison  purposes. 

A  dynamic  simulation  of  a  shipboard  antenna  mounted  on  a  rolling  pitching 
ship  in  a  stop-action  type  of  presentation  is  programmed.   Entering  sea 
state  and  direction  in  the  data  format  causes  the  simulation  to  operate.   Sea 
motion  is  resolved  into  ship  motion  and  ship  motion  into  antenna  parameter 
variation.   Patterns  are  computed  and  displayed  at  10  degree  intervals  of 
wave  period.   Sea  state  0  must  be  entered  to  by-pass  the  dynamic  simulation 
if  it  is  not  desired. 

Appendix  F  is  operating  instructions  for  use  of  the  program  implemented 
at  the  Computer  Graphics  Laboratory,  U.  S.  Naval  Postgraduate  School, 
Monterey,  California.   Figure  A- 16  is  a  schematic  of  the  graphics  computer 
system  at  this  facility. 


2.   PROCESSOR  DESCRIPTION 


A  processor  flow  chart  is  presented  in  figure  A-2.   Processor  operation 
and  interaction  is  described  below. 

A.  The  Parameter  Format  Processor  initializes  the  display  graphics  and 
text  data  blocks  and  displays  the  text  format  for  parameter  and  program 
options  commands  input. 

B.  The  Parameter  and  Options  Input  Processor  is  used  to  enter  problem 
parameters  and  program  option  commands  using  the  format  provided  by  the 
previous  processor.   The  parameter  and  options  input  processor  requires 
entry  of  all  parameters  each  time  the  program  is  initialized.   All  other 
utilizations  of  this  processor  require  changing  only  individual  parameters 
as  desired.   If  the  reinitialization  option  is  selected  by  the  operator, 

the  parameter  format  processor  is  branched  to  from  the  parameter  and  options 
input  processor. 

C.  The  Pattern  Manual  Entry  Processor  displays  a  blank  graphics  data 
block.   By  manually  editing  this  data  block,  the  operator  may  draw  a 
pattern  that  will  be  displayed  with  all  subsequently  computed  patterns. 
Erasing  this  manually  entered  pattern  must  be  done  by  reinitializing  in 
the  parameter  and  options  command  processor.   If  no  manual  pattern  is 
desired,  this  processor  may  be  terminated  without  entry  being  made. 

D.  The  Environmental  Constants  Processor  computes  the  values  of 
problem  constants  that  are  functions  of  antenna  parameters  and  environmental 
conditions  and  not  dependent  upon  observation  angles. 

E.  The  Input  Resistance  Processor  computes  a  value  for  input  resistance 
of  the  antenna  entered  in  the  parameter  and  options  input  processor.   If 

the  equations  in  the  gain  processor  assume  a  nominal  value  of  input  resistance, 
a  value  of  1.0  is  assigned  to  input  resistance. 

F.  The  Observation  Angle  Constants  Processor  computes  values  of  problem 
constants  that  are  functions  of  observation  angles  for  those  observation 
angles  for  which  the  antenna  gain  is  to  be  computed. 

G.  The  Gain  Processor  computes  the  gain  of  the  antenna  selected  in  the 
parameter  and  option  command  input  processor  at  the  selected  zenith  angle 
all  integer  values  of  azimuth  angle  from  1-360  degrees,  and  the  selected 

azimuth  angle  and  all  values  of  zenith  angle  from  1-90  degrees.    These  two 
gain  vectors  are  the  horizontal  and  vertical  gain  patterns. 

H.   The  Normalize  and  Max  Gain  Processor  selects  the  maximum  value  of 
gain  from  both  horizontal  and  vertical  linear  patterns  and  normalizes  both 
patterns  with  respect  to  this  maximum  value.   This  operation  is  required 
to  scale  patterns  for  graphics  display.   The  absolute  value  of  maximum  gain 
is  computed  and  its  log-io  taken.   This  value  is  displayed  in  the  text  data 
format. 

I.   The  Log  Gain  Processor  operates  if  the  operator  manually  selects 
the  log  gain  option  in  the  parameter  and  options  command  input  processor. 
The  horizontal  and  vertical  linear  patterns  are  converted  to  logarithmic, 
patterns  with  a  30  db  range  of  (lOlog^gmax  gain)  to  (lOlog-^gmax  gain)  . 
These  patterns  are  renormalized  by  the  log  gain  processor. 


J.   The  Pattern  Display  Processor  is  a  two  part  processor  which  displays 
the  horizontal  and  vertical  patterns  at  the  graphics  terminal. 

K.   The  Pattern  Save  Processor  is  a  two  part  processor  which  operates 
if  the  horizontal  save  and  vertical  save  option  are  selected  by  the  operator 
in  the  parameter  and  option  command  input  processor.   They  may  be  indepen- 
dently selected.   This  processor  transfers  the  pattern  currently  displayed 
by  the  display  processor  to  storage  in  the  digital  machine  in  a  save  array. 
Entering  a  pattern  in  a  save  array  destroys  the  pattern  previously  saved  so 
care  must  be  exercised  to  bypass  this  processor  if  saving  the  pattern  for 
several  compute  cycles  is  desired. 

L.   The  Display  Saved  Patterns  Processor  operates  when  the  recall  option 
has  been  selected  by  the  operator.   The  processor  recalls  the  patterns  saved 
in  the  save  array  and  displays  them  at  the  graphics  terminal.   Operation  of 
this  pattern  does  not  destroy  data  in  the  save  array. 

M.   Dynamic  Processor.   This  processor  computes  and  displays  a  simula- 
tion of  shipboard  whip,  dipole  and  sloping  longwire  antenna  patterns. 
Entry  of  an  integer  larger  than  0  in  ISEA  will  cause  this  processor  to 
operate.   The  processor  computes  sinusoidal  ocean  waves  with  magnitude  de- 
pendent upon  sea  state.   Ship  roll  and  pitch  which  are  functions  of  ship 
type,  sea  state  and  relative  direction  of  the  seas  are  computed.   Parameter 
variations  caused  by  ship  motion  are  computed  and  the  normal  compute  loop 
entered  with  the  modified  values  of  antenna  parameters.   The  patterns  are 
computed  and  displayed  and  the  gain  at  the  0 '  and  <J>'  inputed  in  the  parameter 
input  processor  is  displayed  under  SIGL  in  the  text  data  format.   The  ocean 
model  is  re-entered.   The  ship  roll  and  pitch  cycles  are  divided  into  36 
discrete  steps  and  a  pattern  computed  and  displayed  for  each  step.   The 
display  will  be  a  stop-action  type  display  of  antenna  pattern  vs.  time. 
Entry  of  0000  under  ISEA  will  cause  this  processor  to  be  bypassed. 


3.   PROCESSOR  FUNCTIONAL  DESCRIPTION 

Figures  A- 3  thru  A-15  are  functional  flow  diagrams  of  processors. 
Equations  for  the  ocean  model,  gain,  and  input  resistance  used  in  the  gain 
and  input  resistance  processors  are  included  in  Appendix  E.   The  sources 
for  gain  and  input  resistance  equations  are  ESSA  Technical  Report  ESSA- 
ERL-110-ITS  78,  A.F.  Barghausen,  J.  W.  Finney,  L.  L.  Proctor,  L.  D.  Schultz, 
May  1969  and  ESSA  Technical  Report  ESSA-ERL-104-ITS  74,  M.T.  Ma,  L.C. 
Walters,  April  1969.   A  listing  of  the  Fortran  Program  used  to  implement 
the  program  is  Appendix  D. 


4.   EXTENSION  OF  THE  PROGRAM 

The  program  may  be  extended  to  compute  patterns  for  other  types  of 
antennas.   Adding  antennas  may  be  accomplished  by  inserting  an  input 
resistance  branch  in  the  input  resistance  processor  and  a  gain  branch  in 
the  gain  processor.   If  additional  parameters  are  required,  the  parameter 
format  must  be  changed  to  accept  them.   Multi-element  antennas  such  as 
Yagi  or  Log  Periodic  will  have  a  mutual  impedance  matrix;  the  terms  of  this 
matrix  may  be  evaluated  using  the  mutual  impedance  equations  in  the  dipole 


branch.   Specific  changes  required  to  add  antennas  to  the  USNPGS  implemen- 
tation are  as  follows: 

1.  Statements  134  and  135  may  be  changed  to  new  parameter  names. 

2.  After  statement  162,  add  DECODE  statements  for  new  parameters. 

3.  New  constants  statements,  if  any,  should  be  inserted  between 
statements  198  and  218. 

4.  After  statement  229  in  the  input  resistance  processor,  insert 
IF(ANTN.EQ.9)  GO  TO  1900. 

5.  After  statement  771  in  the  input  resistance  processor,  add  1900 
INPUT  RESISTANCE  BRANCH  STATEMENTS 


GO  TO  2000 

6.  In  the  gain  processor  after  statement  272,  insert  IF(ANTN.EQ.9) 
GO  TO  900. 

7.  In  the  gain  processor  after  statement  771,  insert  900  ANTENNA 
GAIN  STATEMENTS 


GO  TO  42 

The  dynamic  simulation  is  available  for  whip,  sloping   longwire  and 
vertical  whip  antennas.   Simulation  of  other  antennas  aboard  ship  may  be  made 
by   rewriting  gain  equations  to  allow  arbitrary  orientation  of  the 
antenna.   Orientation  variations  are  available  in  the  ocean  model  and  the 
dynamic  simulation  can  then  be  made. 
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Parameta  Format  Processor 
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Figure  A.  5 
Pattern  Manual  Entry  Processor 
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Low  Gain  Processor 
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Figure  A.  12 
Pattern  Display  Processor 
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Figure  A. 14 
Saved  Pattern  Display  Processor 
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APPENDIX  B 


EXAMPLE  PATTERN  COMPUTATIONS 

Example  pattern  calculations  for  the  seven  antennas  programmed  are 
presented  in  this  section.   Patterns  were  computed  for  typical  parameter 
values  for  each  antenna.   Computation  of  effects  of  parameter  and  environ- 
ment variations  as  well  as  the  use  of  program  control  options  are  demons- 
trated.  The  text  input  required  to  compute  each  pattern  is  presented  with 
a  CRT  photograph  of  the  pattern  computed.   The  USNPGS  user  may  use  the 
text  input  in  conjunction  with  the  user  instructions  of  Appendix  F  to 
learn  program  use. 


Figure  B-21  is  a  film  strip  of  the  36  images  that  comprize  the 
dynamic  simulation  of  a  shipboard  vertical  whip  antenna  in  a  state  5 
sea  from  045  degrees  relative  to  ship's   bow.   Figure  B-22  is  a  dynamic 
simulation  of  a  horizontal  dipole  in  the  same  sea  conditions.   The  images 
of  the  dynamic  simulation  are  computed  at  10  degree  intervals  of  the  ship's 
roll  and  pitch  period.   Figures  B-21  and  B-22  should  be  scanned  down  columns 
and  from  bottom  of  left  column  to  top  of  right  columns. 
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FIGURE  B.l 


ANTN 

0001 

LENG 

01.0 

HGHT 

01.0 

PHIP 

0000 

THEP 

0090 

FREQ 

150. 

EPSL 

10.0 

SGMA 

0.01 

PHI 

0000 

THET 

0080 

FARM 

0000 

_STH 

0000 

ISTV 

0000 

IRCL 

0000 

HGTT 

0000 

ALPH 

0000 

GAIN 

ISEA 

0000 

ICRS 

0000 

SIGL 


Comments:   Manually  entered  A/2 

Vertical  dipole  pattern;  height  A/2 
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FIGURE  B.2 
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Comments : 


a/2   dipole;   A/2    height;  good  ground;  *=  0.   9'=90, 
f=150  mhz,  £=1.0  h=1.0;  overlay  manually  entered  pattern; 

Observation  angles   $=0,  0  =  80 
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FIGURE  B.3 
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Comments :   Erase  manual  trial  pattern 
no  pattern  computed 
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FIGURE  B.A 
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Comments 


A/2   dipole;  X  height;   good  ground;   ^'=0,  0'=90, 
f=  150  mhz,  1=    1.0,  h=2.0;   Observation  angles  <t>  =0 , 
0=80;  save  patterns 
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FIGURE  B.5 
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Comments 


2  A/3   dipole;  A  A/3  height,  good  ground,   <J>'=0,   0'-90, 
f=225  mhz,  £=1.0,  h=2.0;  observation  angles,  4>=0 ,   0=80. 
The  effects  of  changing  frequency  are  shown  here. 
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FIGURE  B.6 
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Comments 


Recall  pattern  4  and  overlay  on  pattern  5.   The  use  of  save 
and  recall  options  are  shown  in  this  example.   The  options 
are  used  to  compare  the  A  dipole  (inside  and  2  lobe  pattern) 
with  4/3X  dipole  (outside  and  3  lobe  pattern). 


FIGURE  B.7 
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Comments:  Two  thirds  wave  length  dipole;  Four  thirds  wave  length 
height;  good  ground  ,4> '  =0  ,0  '=90  ,  f=225  mhz,  £  =  1.0,  h=20; 
observation  angles   4>=0   0=80  ;Log  patterns  30  db  scale. 

Log  pattern  option  is  used  to  study  side  lobe  structure. 


34 


FIGURE  B.8 
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Comments:   Two  thirds  wavelength  dipole;  four  thirds  wave  length  height, 
good  ground,  <j>'=0,   0'*45°  (tiltangle),  f=225  mhz,   £=1.0 
h=2.0,  observation  angles  $=Q      0=80;   Tilted  dirole.   The 
effect  of  tilt  on  dipole  radiation  patterns  is  demonstrated  here, 
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FIGURE  B.9 
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Comments:   Quarter  wavelength  whip;  good  ground,   f=30  mhz,  £=25 
observation  angles   (J)=0 ,  0=80 
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'1GI  RE    3.  10 


ANTN 

0001 

LENG 

02.5 

HGHT 

00.0 

PHIP 

0000 

THEP 

0000 

FREQ 

030. 

EPSL 

04.0 

SGMA 

0.01 

PHI 

0000 

THET 

0080 

PARM 

0000 

ISTH 

0000 

ISTV 

0000 

IRCL 

0000 

HGTT 

01.0 

ALPH 

0000 

GAIN 

ISEA 

0000 

ICRS 

0000 

SIGL 


Comments:   Quarter  wavelength  whip;  poor  ground,  f=30  mhz,  £=  2.5, 
observation  angles  *=0 ,  0=80.   The  effects  of  changes  in 
reflecting  "round  are  shown  in  this  example.   The  ground 

change  from  good  ground  to  poor  ground  causes  a  decrease  in 
gain  of  2db  and  a  slight  increase  in  0  of  max  .radiation. 
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FIGURE    B.ll 
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Comments :       A/4    vertical    whip   with 
half   wave    length  ground    radial    screen 
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FIGURE    B.12 
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Comments :       A/4    vertical   whip   with   wave    length 
ground    screen.       Increasing    screen   narrows    and 
depresses    vertical    pattern. 
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FIGURE  B.13 
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Comments:    Inverted  L,  horizontal  run  two  wavelengths,  vertical  run 
one  wavelength;  good  ground;  h=10.0,   £=20.0,  f=30  mhz; 
observation  angles   $=60,  0=80. 
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FIGURE  B.14 
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Comments 


Inverted  L,  horizontal  run  two  wavelengths,  vertical  run 
one-half  wave  lenth;  good  ground;   h=5.0,  £=20.0,  f=30  mhz; 
observation  angles  0=90,  0=80.   The  effects  of  change  in 

vertical  run  length  are  shown  here. 
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FIGURE  B.15 
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Comments:   Sloping  Long  wire;  two  wavelengths;  good  ground;   4>=0 , 
0'=45°,  £=20,   f  =  30  mhz,   observation  angles  <f>=0 ,  0=80, 
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FIGURE  B. 16 
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Comments:   Sloping  Longwire;  two  wavelengths;  good  ground;  $'=0,0'  =60, 
£=20.0,  f=30  mhz;  observation  angles  $=0 ,  0=80.   This  set  of 
two  examples  demonstrates  the  effect  of  variation  of  tilt 
angle  on  radiation  patterns. 
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FIGURE  B.17 
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Comments:  Horizontal  rhombic;  three  wavelength  sides;  good  ground; 
one  wavelength  height;  £=30.0,  h=10.0,  a=30° ,  f=30  mhz; 
observation  angles   $=0,  0=80. 
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FIGURE  B.18 
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Comments:   Horizontal  rhombic,  three  wavelength  sides;  good  ground, 

one  wavelength  height;      £=30.0,      h=10.0      a=45°,      f=30mhz; 
observation    angles   $=0,    0=80.      These   last    two    computations 
show   clearly   how   the   program  may   be      used    to   synthesize 
antenna   systems.      A  non-optimum  a   is    compared    to    the  optimum 
for   a   given  h ,£   etc.      Since    this    antenna   is    fairly    difficult 
to   build,    the   use   of    the  program   to   synthesize    the   optimum 
is   well   justified. 
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FIGURE  B.19 
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Comments:  Vertical  half  rhombic,  three  wavelength  sides,  good  ground, 
£=30.0,  a=30°,  f=30mhz;  observation  angles  $=0,  9=80.  This 
antenna  has  a  major  lobe  at  approximately  40°  elevation  and 
may  be  suited  for  propagation  conditions  requiring  high 

elevation  lobes. 


FIGURE  B.20 
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Comments:    Vertical  half  rhombic,  three  wavelength  aides,  good  ground, 
£,=30.0,  a=45°  ,   f  =  30  mhz;  observation  angles  (j)=0 ,  0=80. 
Increasing  a  splits  the  energy  into  a   high  and  low  lobe  and 
decreases  gain  from  the  previous  case. 
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Vertical 


Whip  Pattern 


Small  Ship 

Sea  State  5  from  045°R 


Figure  B-21 
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Figure   B-21 
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Figure   B-21 
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Vertical  Whip 
Horizontal  Pattern 


Figure  B-21 
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Figure   B-21 
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Figure  B-21 
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Half  Wave  Dipole 
Vertical  Pattern 


Small  Ship 

Sea  State  5  from  0A5°R 


Figure  B-22 
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Figure  B-22 
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Figure   B-22 
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Figure  B-22 
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Half  Wave  Dipole 
Horizontal  Pattern 


Figure  B-22 
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Figure  B-22 
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Figure   B-22 
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APPENDIX  C 

ANTENNA  GEOMETRY  AND  GAIN  AND  INPUT  RESISTANCE  EQUATIONS 

This  appendix  details  the  antenna  geometry  of  the  programmed  antennas 
and  the  equations  used  in  the  gain  processor  and  input  resistance  processor. 
The  sources  of  equations  and  geometry  are  ESSA  Technical  Report  ESSA-ERL- 
110-ITS  78  and  ESSA  Technical  Report  ESSA-ERL-104-ITS  74.   The  spherical 
coordinate  system  used  ro  describe  antenna  patterns  is  the  IEEE  standard 
and  is  shown  in  figure  C-l   Antenna  geometry  for  the  antennas  programmed 
is  shown  in  figures  C-2  thru  C-9 . 

The  definitions  of  the  terms  used  in  antenna  equations  are  as  follows: 

I   =   length  of  a  unit  radiator  in  meters 

h  =  height  of  antenna  feed  point  above  ground  plane 

0""  Tilt  angle  of  the  antenna  axis  measured  from  the  zenith 

A "=   Tilt  angle  of  the  antenna  axis  measured  from  the  horizontal 

0  =  Observation  zenith  angle 

A  =  Observation  elevation  angle 

a  =   Apex  half  angle 

a^-   Complement  of  apex  half  angle 

<|>  =  Observation  azimuth  angle 

£r=  Dielectric  constant  of  ground  plane 

o  =  Conductivity  of  ground  plane 

R,  =  Complex  ground  reflection  factor  for  a  horizontally  polarized  wave 
h 

CH  =  Magnitude  of  horizontal  reflection  factor 

4V-=  Phase  of  horizontal  reflection  factor 

R  =  Complex  ground  reflection  factor  for  a  vertically  polarized  wave 

CV=  Magnitude  of  vertical  reflection  factor 

T  =  Phase  of  vertical  reflection  factor 
v 

R  ,,=  Complex  horizontal  reflection  factor  evaluated  for  normal  incidence 
h 

R  ^=  Complex  vertical  reflection  factor  evaluated  for  normal  incidence 
v-" 

f  =  Frequency  in  mhz 
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A   =  wave  length  in  meters 

a  =   radius  of  ground  screen 

c  ■   radius  of  ground  screen  wire 


Equations  for  quantities  that  are  common  to  all  antennas  programmed 
are  as  follows: 

X   ,=      3.0   x  108 


k     =     2II/A 
k2»    1  I 


r    =  cose  -  k 

*2 


L      -  j      1.8  x   104g    j  * 

[-ft-)7 


cose  +  k    fi  -A    siNe  \      | 

k2L  vk2     )  J 

r^  =  cose  -  k2  |"i  -/k      siNe  \2|  1/2 


cose  +  k    [i  -/k      sino  V  1 1/2 

IT  L   \^2       /  J 


Rv'  =  k2  -  k  V  =  k  -  k2 

k2  +  k  k  +  k2 


51  =   COS  (ik  -  2Kh  SINA  ) 

n 

52  =   SIN  (V,  -  2Kh  SINA  ) 

h 

53  =   COS  (Y   -  2Kh  SINA  ) 

v 

54  =   SIN  (T  -  2Kh  SINA  ) 

v 
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C-  1  ARBITRARILY  TILTED  DIPOLE 

The  equations  presented  are  for  a  thin,  single  element,  center  fed 
dipole  arbitrarily  oriented  above  a  flat  ground  plane. 

GI  =   COS  (l/2k«,  (SINA  SINA'  +  COSA  COSA '  SIN  4>  )  -  COS  (l/2k&) 
1.0  -  (SINA  SINA'  +  COSA  COSA'  SIN<£)2 

DI  =   COS  (l/2k&  (COSA  COSA'  SIN  <j>  -  SINA  SINA'  )-  COS  (l/2kfc  ) 
1.0  -  (COSA  COSA'  SIN<j>-  SINA  SINA'  )2 

E^   =  (COSA'  SIN  4>  SINA  -  SINA'  COSA  ) -GI 

-  (COSA'  SIN  <J>  SINA  +  SINA'  COSA  )-DI-CV'S3 

E^.   =  COSA'  COS  <}>  (GI  +  DI-CH-S1) 

E    =  (COSA'  SIN  $  SINA  +  SINA'  COSA  )-DI'CV»S4 
J2 

E 
l        =   COSA'   COS  (f>  •DI-CH'S2 
*2 

Gain  =  120-  (E   2  +  E   2  +  E   2  +  E   2)  /Rin 


S  •  2h 


S   =  S  SING'  COS  d)1 
x 

S   -  S  SING'  SIN  d>' 

y 

S   -  S  COSQ 
z 


p.[ 


S  2  +  (Y  +  S  )2    1/2 
x       o     y    J 


Y   =2  COSA '  h 
o  — 

A   \    for  mutual  impedance 

Z   -  2  SINA'  h 
I 


for  self   impedance 


Y      -       /H x   10    3    •     n/X 
o 

Z      =  .0 
o 

r     =[p2+(Z     +S)2]1/2 
I  o  z      J 

r     =      p2  +    (Z     +  S     +  1/2*   )2        1/2 
1        L  o  z  J 

r,   =      p2  +    (Z     +  S      -   1/2*    )2  J     1/2 
Z        ^  o  z  J 
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SR  =    1/r   SIN    (2nr) 


SRI  =   l/r1   SIN    (2nr1) 


SR2   =    l/r2   SIN    (2nr2) 


FACR  =    2*SR  COS    (IU) 
CR  =    1/r  COS    (2IIr) 


CR1   =   l/r1   COS    (2nr1) 


CR2   =   l/r2   COS    (2nr2) 


FACX  =    2«CR«COS    (IU) 


Z. .    =    (R. .   +   j    X. .    ) 
iJ  ij  ij 


h,-»fl'2l\l( 


SR1.(S      +   Z      +   £)+   SR.    (S      +  Z      -    £  ) 
z  °T  2z  o        y  ' 


-   FACR-(S     +  Z   )  )'(S    : 

Z               O      /           X 

rsiN2n(£  -  Isl )  1 

2                      dS 

L       s         J 

r  ll1  (r      / 

X.  .= 

-30 

1      k  CRWv 

■] 


2+YS      +   S    2)      +   S      (FACR  -   SRI   -    SR2) 
o   y  y  z 


}' 


-  A    )    -   FACR    (S      +  Z    )  )• 
2 


:') 


] 


(S2  +  YS      +  S    2)  I  +   S      «(FACX  -   CR1-CR2) 
x  o  y  y  z 


! 


SIN2II    (£   -    |S    ) 


] 


dS 


Z..,  =  self  impedance 
Z   =   mutual  impedance 


R.   =  R 
in    11 


+  Real 


Z21  (Rh'  C0SA'   +  J  R  '  SINA?)  (COSA'  -  j  SINA 


"] 
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C.2   VERTICAL  MONOPOLE 

Vertical  whip  antenna  equations  are  for  a  base  loaded  vertical  whip 
above  a  flat  ground  plane. 

S3  =  COS  (4>  ) 

s4  -  sin  (y  ) 

v 
A  =  COS  (k£  SINA)  -  COS  (k£) 
B  =  SIN  (kl   SINA)  -  SINA  SIN  (k£) 

Gain  =  )  [a»  (1+CV-S3)  +  B-CV'S4]2  -I-  [  A-CV  S4  +  B«  (1-CV*S3)  ]  2 

'  R.   COS  A 

in 

R   -  15  <[2  +  2  COS  (2k£)J»[en  (2kft)  +  Y  -  Ci  (2k£)]  - 

COS  (2kA)-[jln  (4k£)  +  Y-  Ci  (4kfc)  ]  - 

2  SIN  (2k£)'Ptf  +  Si  (2kA)  1  +  SIN  (2k£)«U_+  Si  (4k£)   I 

I2  J  I1  J) 

Y  =  .577 


Ci(x)  =  - 


cost 


dt 


Si(x)  =  - 


SINt 


dt 


'x 
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C.3  VERTICAL  WHIP  WITH  GROUND  SCREEN 

The  vertical  whip  with  ground  screen  is  a  single  monopole  above  a  flat 
ground  with  a  radial  conductor  ground  system  consisting  of  N  equally  spaced 
radial  conductors.   A  value  of  120  is  used  for  N  and  1  cw  wire  is  assumed. 

53  =  COS  ^ 

v 

54  =  SIN  i|j 


A  =  COS  (k£SINA)  -  COS  (k£) 


B  =  SIN  (k£SINA)  -  SINASIN  (kl) 


Gain  = 


2    2 
[A-(l  +  CV»S3)  +  B-CV-S4]2  +  [A-CV-S4  +  B- (1-CV-S3) ]2  ^A3  +  B3j 

R. 

in 


A3+JB3  - 


n  = 


rka  2   2  2.1/2 

1  -  nSINS     [  e  nX  }  -  e  JXC0S(k£)]  J1(xSIN9)dx 

120n  SIN(k£)  [COSCkJICOSe)  -  C0S(k£)] 
1/2 


f       j^      1 

o  +  j  a)  e  J 


R,    =  R.      of  vertical  whip    (C.2) 
1  in 


«  - 


4ns IN    Qui) 


12kA 


-J2WW4 


e^^Ei  ["J2k(r0  +   £)  ]     +  e  Ei  [-J2k(r0"£)  ^ 


+   2C0S    (k£)    Ei[-j2ka]   +  4   C0S(k£)    Ei  [-jkr^ 


-   4C0S(k£)e  JkS,EI 


C0S(k£)ejk£E: 


AZ, 


■n[s-J 


L  [-JkC^  -   H)l     -  4 

r/  2        2    1^2 

I    -jk(p  +  «, 1        _e"jkp    C0S(k£)l 


i   r-jk(r1+£)J  | 


2IIp    SIN    Qui) 


dp 


Ei   (+  jx)    =  Ci(x)    +  j    Si(x) 

,    2^fl2,l/2 
r     =    (a  +1   ) 
o 


rx   =   a  +    (a2+£2)1/2 


CT+]W£j 

n       fj240n2pUn^  \ 
'e     \         NA  /       VNc   / 


R.      =   R,    +  Real    (AZ,    +  AZ.) 
in  1  12 
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C.  4   INVERTED  L 


The  inverted  L  antenna  equations  are  for  a  long  wire  antenna  that  is 
base  loaded  and  arranged  in  an  inverted  L  configuration. 

A  -  COS(k£)  COS  (khSINA)  -  SINA  SIN  Qui)    SIN  CkhSINA) 

-  COS  (k  (h+fc)) 

B  =  SINA  SIN  (k£)  COS  (kh  SINA)  +  COS  (k£)  SIN  (kh  SINA) 

-  SINA  SIN  (k  (h+O  ) 

GI  =  SIN  (kS.  COSA  SIN  <})  -  COSA  COS  $  SIN  Qui) 
GR  =  COS  (k£  COSA  SIN  <J)  -   COS  Qui) 

-S3)  +  GI-CV-S4)  g 


0 


[SIN  4>  SINA  LgR(1, 
1.0  -  C0S2A 


SIN2  <j> 


[ 


A  (1.0  +  CV«  COS  i  +  B  CV.SINi^ 
_v v 

COSA 


] 


SINj  SINA  (GI  (1.0  -  CV'S3)  -  GR'CV»S4 
1.0  -  C0S2A  SIN2  ij> 


B  (1.0  -  CV'COS  Jb  )  +  A'CV'SIN  ^ 


COSA 


♦' 


cos  4 


1  -  cos2a  COS  4> 


] 

,|lGR(l.  0  +  CH» 


SI)  -  GI-CH-S2 ] 


+  [GI  (  1.0  +  CH-S1)  +  GR-CH-S2  ] 
Gain  =  30-0  [|e/+  |e/  ] 


) 


Rin 


Si(x)  =  I 
x 


dx 


Rin  =  60 


L.41  +  I 


n  I  21  \   +  SIN  (2kg,) 
\  X  /      2k£     J 

(Jin  (2kh)  +  1.270  -  Ci  (4kh) 

In   (2kh)  +  0.577  -  Ci  (2kh)) 

1/2  Si  (4kh)  -  Si  (4k£)  I 


+  30.0  I-   1/2  COS  (2kh)|Un  (2kh)  +  1.270  -  Ci  (4kh) 
+  (1.0  +  COS  (2kh))|  ^n  (2kh)  +  0.577  -  Ci  (2kh)) 
-  SIN  (2kh)   1/2  Si  (4kh)  -  Si  (4k£) 
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C.5   SLOPING  LONG  WIRE 

The  equations  for  this  antenna  are  for  a  base  loaded  longwire  antenna 
arranged  in  a  sloping  configuration;  the  antenna  zenith  angle  may  be  assigned 
values  of  0  thru  90  degrees. 

t  i 

CIG  =    COS  [k£  (SINA  SINA   +  COSA  COSA   COS*)]  -  COS  (ki) 

i  i       o 

1.0  -  (SINA  SINA   +  COSA  COSA   COS  *) 

i  » 

SIG  =    SIN  [kg.  (SINA  SINA   +  COSA  COSA   COS  *)  ]    -(SINASINA  '+COSAC0SA  'COS*)  SIN(k&' 

i  i       2 

1.0  -  (SINA  SINA   +  COSA  COSA   COS*) 


CIGP  =    COS  [kfc)  (  COSA  COSA   COS*-  SINA  SINA   )]-  COS  Qui) 

i  i  2 

1.0  -  (COSA  COSA   COS*-  SINA  SINA  ) 


1  2 
1-0  -  (COSA  COSA   COS*-  S-INA  SINA  ) 

f  i 

+  (SINA  SINA   -  COSA  COSA   COS  *)  SIN  (ki) 


r 

SIGP  =    SIN['k&)  (  COSA  COSA   COS*-  SINA  SINA  )] 
L  LO  -  (COSA  COSA   COS  *  -  J 

]      1 

E    -   -COSA   SIN  *  J  CIG  +  CH  '(CIGP  COS  *   -  SIGP  "SIN  it  ) 
E      =  -  COSA   SIN*   SIG  +  CH  (CIGF'SIN*   SIGP  COS*) 

!  I  » 

E     =  CIG  (COSA   COS  *  SINA  -  SINA   SINA)  +  CV  (COSA   COS  *  SINA 
1 

+  SINA   COSA)*  I  CIGP  -COS  tL  -  SIGP*  SIN  $   1 

F  '  '  *  ' 

02    =  SIG  (COSA   COS  *  SINA  -  SINA   COSA)  -  CV  (COSA   COS  *  SINA+  SINA   COSA) 
IciGP'SIN*  +  SIGP'COS  * 

Gain     =    30  |E.    2  +   E  ,     2     +  E      2     +  E        2  1/    Rin 

L  *1  *2  *i  *2     J 

Rin  =   30  { 1/2 [an  (ki)   +    .577  +  Ci    (4k£)  | 
+  0,693  +  COS    (kO     [COS   k£    Qln  (k£)   +.577 

-  2.0   Ci    (2k£)    +  Ci    (4k£))    -   SIN    (k£)    (Si    (4kjO 

-  2.0   Si    (2kH)  I   ( 

/oo       J       ) 
x^-  dt 

"J!  2?kat 


Si(x)    = 
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C.  6   TERMINATED  SLOPING  VEE 

The  terminated  sloping  vee  equations  are  for  two  sloping  longwire 
antennas  arranged  in  a  vee  configuration.   The  feed  point  is  the  apex  of  the 
vee.   The  elements  are  fed  180  degrees  out  of  phase.   The  elements  of  the  vee 
are  terminated  in  370  ohm  non  inductive  resistors. 


COS  t|) 


COS  \\i 


COS  i- 


COS  i< 


cos  i 


cos  i 


cos  i 


cos  i 


SINA  SINA 

SINA  SINA 

-SINA  SINA 

-SINA  SINA 

COSA  SINA 

COSA  SINA 

-COSA  SINA 

-COSA  SINA 


+  COSA  COSA 

+  COSA  COSA 

+  COSA  COSA 

+  COSA  COSA 

+  SINA  COSA 

+  SINA  COSA 

+  SINA  COSA 

+  SINA  COSA 


COS 
COS 
COS 
COS 
COS 
COS 
COS 
COS 


4»  -  <* ) 

<()+  a  ) 

<J>  -  a  ) 

4>+  o  ) 

t»-  a) 

cf  +  or  ) 

♦  +  -a) 


U  =  k£  (1.0  -  COS  ij,  J 
i  i 


i  =  1,  2,  3,  4 


A  =    COS ^  _  (COS  (U  )  -  1.0)     -   COS i   .  (COS  (U0)  -  1) 


+  cv . 


COS  i ■        S3(C0S  (U4)  -  1-0)  S4  +  SIN  (U^) 


-  COS  ip      (COS  (U3)  -  1-0)  S3  +  SIN  (IL)  S4] 


B  = 


COS  ty     SIN  (IL)   -  C0Si|;  SIN  C^) 


+   CV    J  COS  i>  c    I  SIN    (U3)    S3   -    (COS    (U3)    -   1.0)    S4    j 


( 


+   COS  i  ,       I    (COS    (U,)    -   1.0)    S4      -   SIN    (U)    S3 
oi4  4 


) 


+   C0Si|<,       (    (COS    (IL)    -    1.0)    S4      -  SIN    (U)    S3 
6      V  4 4 

x  U 
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C  =        SIN    U  +  a)     (COS    (U2)    -    1.0)      -   SIN    (  cj>  -  a )  CCOS    Cljj)    -   1.0  ) 


,A)    -   1.0)      -   SIN    C4>  -a)    (COS    CU3) 


+   CH    J   SI    ("siN    (<fc+a)    (COS    0} 

(     L  ^ 

-S2[siN    (<|>-a)    SIN   (U  )  -      SIN    (  <fr  +  a  )    SIN  (1^)1 


U., 


1.0)    j 


D      =        SIN    U-a)SIN    (U   )         -      SIN    U+a)    SIN    QJ   ) 


+   CH 


|~SIN    (  4>  -  a  ) 


SIN    (U   )         -   SIN    (<}>  +  a)    SIN    (U   )|.sl 


} 


+        SIN    (  $ +  a 


D 

'3  4 

)    (COS    (U   )    -10)      -     SIN    (<(>-a)    (COS    (U3)    -   1.0   ) 


S2 


Gain  =   0.05 


o  0  9     '  7  2 

A     +  B      +   COSZA  (C     +  D   ) 
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C.  7   HORIZONTAL  RHOMBIC 

The  horizontal  rhombic  antenna  equations  were  developed  under  the 
assumption  of  uniform  current  distribution  of  the  effective  value  of  current, 
The  antenna  is  loaded  at  the  apex  and  terminated  in  dissipation  lines  at  the 
opposite  corner. 


u\    =   1 .0   -  COSA    CSIN  a         COS  <J>  +  COS  a         SIN  cj>  ) 
1  c  c 


U„  =   1  ;0   -   COSA    (SIN  a         COS  $  -   COS  a         SIN  <j>  ) 
2  c  c 


Gain  =  2  .16 


COS  a     SIN  (1/2  kAUjSIN  (1/2  k£U.) 
c ; 1_ 2 

U1U2 


(COS  <j>  -   SIN  o    COS      )2    (CH2   +   1  .0  +   2  .0    -CH-    SI)    + 
SIN2A      SIN2  4>     (CV2   +   1.0-2.0    -CV-    S3) J 
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C.  8   VERTICAL  HALF  RHOMBIC 

The  vertical  half  rhombic  are  for  a  base  loaded  longwire  arranged  in  the 
vertical  half  rhombic  configuration  and  terminated  in  a  40Q-50Q  ohm  non-induc- 
tive resistor. 


SI  =  SIN  (K  (1.0  -  COS  i     )) 


CI  =  COS  (K  (1.0  -  COS  t|<  )) 


S2  =  SIN  (K  (1.0  -  COS  ty     )) 


C2  =  COS  (K  (1.0  -  COS  i|;  )) 


COS  i       =  COSA  COSA   COS  <j>  -  SINA  SINA 


COS  Tj;   =  COSA  COSA   COS  <j>  +  SINA  SINA 


Rl  =  1  -  CI  II  =  SI 


l.o  -  cos  i>  i.o  -  cos  ^ 

R2  =  CI  (1.0  -  C2)  +  S1-S2  12  =  C1-S2  -  SI  (1  .0  -  C2) 

i.o  -  cos  i  1.0  -  cos  i 

R3  =  (1.0  -  CI)  COS  (2k 2,  SINA  SINA)  +  SI  SIN  (2k£  SINA   SINA) 

i  t 

13  =  SI  COS  (2k£  SINA   SINA  )  -  Q-  .0  -  CI)  SIN  (  4IU  SINA   SINA) 


Fl  =  I3-C1  -  R3-S1 
1.0  -  COS  \\) 


F2  =  R3-C1  +  I3-S1 

F3  =  1.0  -  C2 

1.0  -  COS  ty  2 

F4  =     S2 


1  -  COS  4' 
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RB  =  Rl  +  R2  -  CV  [(F2  +  F3)  S3  -  (Fl  +  F4)  S4  ] 

Bl  =  II  +  12  -  CV  [(F2  +  F3)  S4  +  (Fl  +  F4)  S3  ] 

RC  -  R2  -  Rl  +  CV  [(F2  -  F3)  S3  -  (Fl  -  F4)  S4  ] 

CC  =  12  -  II  +  CV  [(F2  -  F3)  S4  +  (Fl  -  F4)  S3  ] 

RA  =  Rl  +  R2  +  CH  [(F2  +  F3)  SI  -  (Fl  +  F4)  S2  ] 

Al  =  II  +  12  +  CH  [(F2  +  F3)  S2  +  (Fl  +  F4)  Si] 

r  '  *2 

Gain  =  0.1    (RB  COSA   COS  <£  SINA  +  RC  SINA   COSA) 

'  '       2 

+  (Bl  COSA   COS  $  SINA  +  CC  SINA   COSA) 

+  (RA  COSA   SIN  <j>  )  2  +  (AI  COSA   SIN  $  )  2 
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Figure  C-l 
Spherical  Coordinate  System 
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L 
H 


9 


ength  of  dipole 

height  measured  al  center 

of  dipole 

Tilt  angle  (measured  in  zenith) 

Tilt  angle  (measured  in  elevation 
above  horizontal) 


<f>      Train  angle  (measured  in  azimuth) 


Figure  C  -2 
Arbitrarily  Tilted  Dipole  Geometry 
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' 


L  Length 


• 


Figure  C-3 
Vertical  Whip  Geometry 
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L       length  of  whip 

H       radius  of  ground  screen  radial  elements 

Figure  C  -  4 
Vertical  Whip  with  Ground  Screen  Geometry 


76 


►    X 


H     -      Height  of  vertical  run 
L     -      Length  of  horizontal  run 


Figure  C-5 
Inverted  L  Geometry 
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L    —    Length  of  Long  wire 

9   —    Tilt  angle  (measured  in  zenith) 

a   —     Elevation  angle 


Figure  C-6 
Sloping  Long -Wire  Geometry 
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370  X2 


L 

a 
9 

H 


Length  of  sloping  element 
Half  element  separation 

Slope  angle  (measured  from  zenith) 
Height  of  load  point 


Figure  C-7 
Terminated  Sloping  Vee  Geometry 
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L   - 

Length  of  side 

H  - 

Height 

a  - 

■    Half  separation 

angle 

q^  load  point 

a  - 
c 

Complement  of 

a 

Figure 

C-8 

Rhombic 

Geometry 
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z 

A 


t 

H 


k 


>x 


L    - 
H    - 

8'- 

a  - 


Length  of  side 
Height 

Tilt  angle  of  side   (measured  from  zenith) 

Half  separation  angle  of  rhombus  at  load 
point 


Figure  C-9 
Vertical  Half- Rhombic  Geometry 
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APPENDIX  D 


PROGRAM  LISTING 

This  appendix  is  a  listing  of  the  Fortran  implementation  of  the  antenna 
pattern  graphics  program.   The  listing  is  preceded  by  a  partial  listing  of 
definitions  of  the  computer  variables  used.   Variables  used  in  graphics 
display  processors  only  were  not  included  in  the  below  list. 


ANTN 

L 

H 

PHIPR 

THEPR 

F 

EPSLN 

SIGMA 

M 

KAY 

PAR 

ISTRH 

ISTRV 

ALPH 

ALP  CM 
DLPRI 

LMDA 

K 

C2 

RHPRI 

RVPRI 

S 


Antenna  Type 

length 

height 

? 
0 

f  in  mhz 

e 
r 

a 
0  in  degrees  (observation  zenith) 

<f>  in  degrees  (observation  azimuth) 
Reinitialize  and  Log  Gain  option  command 
Save  horizontal  pattern  option  command 
Save  vertical  pattern  option  command 


a 
< 

A' 

A 

k 

k2 

V 


R  , 
v 
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sx 

S 

X 

SY 

S 

y 

SZ 

s 

z 

YO 

Y 
o 

ZO 

Z 
o 

ROW 

P 

R 

r 

Rl 

rl 

R2 

r2 

PI 

n 

RIN 

R. 
m 

THETA 

0 

PHI 

* 

KCOS 

COS 

(0) 

RV 

R 

V 

RH 

\ 

SIGHV 

K 

SIGHH 

H 

DELTA 

A 

COSDL 

COS 

(  A  ) 

SINDL 

SIN 

(  A  ) 

SINDP 

SIN 

(a') 

COSDP 

COS 

(a') 

SINPI 

SIN 

(<}>  ) 

COSPI 

COS 

(*  ) 
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ETHT1 

\ 

ETHT2 

\ 

EPHI1 

Eh 

EPHI2 

\ 

) 

G 

Gain  (relative  pattern) 

GAIN 

Max 

Gain 

NORM 

Max 

value  of  linear  gain 

EPHI 

|E* 

,2 

ETHET 

|E* 

I2 

KOS1 

COS 

♦  i 

KOS2 

COS 

*2 

KOS3 

COS 

*3 

KOS4 

COS  i\> 

KOS5 

COS 

^5 

KOS6 

cos  i  6 

KOS7 

COS  t\>  7 

KOS8 

COS  \\> 

COSU1 

COS 

(Ul) 

COSU2 

COS 

(U2) 

COSU3 

COS 

(U3) 

COSU4 

cos 

(U4) 

SINU1 

SIN 

(Ul) 

SINU2 

SIN 

(U2) 

SINU3 

SIN 

(U3) 
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SINU4 
SINAC 

COSAC 

IRCAL 

ISEA 

ICRS 

SIGL 

ADA 

DPHIP 

WAVE 

DLTI 

DLT2 

DLT3 

SIND3 

SINA 

COSA 

VAR 

Z 

RGRAL 

XGRAL 

CEE 

SRFAC 

DLTZ1 

DLTZ2 

CV 
CH 


SIN  (U4) 


SIN  (a  ) 
c 


COS  Co) 
c 


Recall  saved  pattern  option  command 
Sea  State 
Sea  direction 


Log1Q  (G  CM,  KAY)  ) 


n 

A  <j> 
wave 

A 

ei 

A 
0. 


SIN 

(A03) 

SIN 

CA<J>) 

COS 

(Ac))) 

U)t     I 

[wave) 

z 

Ri  -i 

K1J 

Xij 

C 

CA3 

+   JB3) 

AZ. 


AZ 


2 
Rv| 

Rhl 
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Computer  variablesv that  are  identical  to  the  terms  in  Appendix  C  they 
represent,  are  not  listed  here. 

The  following  sub-programs  are  included  in  the  program: 

1.   SUBROUTINE   SINUS  (X,  SC) 

COMPUTER   Si  (x)   =  -      SIN  t 


/SIN 


dt 


2.       SUBROUTINE      KOSINUS    (X,    CC) 

/oo 
COS    t 
x  - 


dt 


3.  FUNCTION  CINC  (x) 

COMPUTER  COS  x 
x 

4.  FUNCTION  SINC  (x) 

COMPUTER  SIN  X 
x 

5.  FUNCTION  AKEX  (x) 

COMPUTER  Ei  (±.jx) 

6.  FUNCTION   ADAE    (x) 

COMPUTER     n 
e 

7.  FUNCTION   ZGRAL    (x) 

COMPUTER  integrand  for  AZ2  =     (   )  dr 

8.  FUNCTION  RESIST  (s) 

COMPUTER  INTEGRAND  for  Rij  =   I  (   )  ds 

9.  FUNCTION  REACT  (s) 

COMPUTER  INTEGRAND  for  Xij  =     (   )  ds 
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APPENDIX  E 


SHIPBOARD  ANTENNA  DYNAMIC  SIMULATION  EQUATIONS 

This  appendix  presents  the  development  of  ship  motion  equations  as 
functions  of  sea  state  and  relative  direction  of  the  sea.   Ship  motion  is 
resolved  into  parameter  variation.   The  values  for  time  varying  parameters 
are  used  in  the  compute  loop  for  the  dynamic  simulation. 

The  ship-ocean  combination  is  modeled  as  follows: 

1.  The  ship  will  roll  sinusoidally  8  degrees  per  sea  state  if  the  sea 
is  on  the  beam,  ie .   from  090   R  or  270   R. 

2.  The  ship  will  pitch  sinusoidally  2.4  degrees  per  sea  state  if  the  sun  is  on 
the  bow  or  stern,  ie.   from  000   R  or  180   R.   (This  represents  a  small  naval 
combatant   ship) 

3.  Sea  state  and  direction  is  resolved  into  ship  motion: 

Vertical  Whip  or  Sloping  Long  Wire 
AA  =  2«L*  SIN  (AG1/2) 
BB  =  2-L-  SIN  (A91/2) 
CC  =   AA2  +  BB2 

DD  =   L2  -  (CC/2)2 

|  A   |  =  2. tan  _1   (CC/2) 
03  DD 

A   -  |A   | 

o3         e  3  ^ 

o'  (t)  =   e'  (0)  -  A      (TILT) 

°3 

SIN  A   =   SIN  AG  SIN  A 

e3  o3 

2    1/2 
COS  AG  =   (1  -  SIN  AG) 

I A  $  (   =  tan"1    SIN  A  fycOS  A  $        (azimuth  change) 
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Wave 

ROU  : 

A0I 

Pitch 

A02 

=  8-SEA-SIN  (wt) 

=  Wave  •  SIN  (CRS.  ^q) 

=  0.3- wave- cos  (CRS.  %<,) 

Pitch 


Dipole 


B\t)  =  9l  (0)  -   l9z 


(tilt) 


H   (r)    =H    (O)-cos   (A0l)-cos    (A0 1  )    (height) 
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if:      A0  >  0  and  AQ  >   Q                                     A  $  =    |A$| 
1                                 2 

A0  <  Q  and  AQ  >   0                                       A'$=-|Acf| 

1  2 

AQ  >  0  and  A0  <    0                                       A  $  =    (II- |  A  <p  \) 

A0  <  0  and  A0  <,   0                                       A  $  =   —  Cn   -|A<j>    |) 

*     (t)    =  <J>'    CO)  -   A  i 


Definition  of  terms: 
A    -   ship  roll* 

ei 

A    -   ship  pitch 

S2 
A   -   tilt  of  whip  or  long  wire  caused  by  ship  motion  (whip) 

G3 
wave  -  sinusoidal  wave 

to    -  wave  radar  frequency 

0  (t)-  antenna  tilt   (Dipole) 

h  (t)-  antenna  height  (Dipole) 

A  4>   -  variation  in  antenna  train  caused  by  ship  motion  (whip) 

<(>  (t)  -  antenna  train  (whip) 

CRS  -  Direction  of  sea  relative  to  ship's  bow 


121 


APPENDIX  F 


OPERATING  INSTRUCTIONS  FOR  U.S.  NAVAL  POSTGRADUATE  SCHOOL  GRAPHICS  COMPUTER  LAB 

This  appendix  gives  step  by  step  operating  instructions  required  to  use 
the  antenna  patterns  graphics  program  at  the  Naval  Postgraduate  School.   Use 
of  the  graphics  library  program  "GATED"  and  computer  light-off  procedures  are 
covered  in  the  operators  manual  and  laboratory  memoranda  and  are  not  included 
in  this  appendix. 

1.  Light-off  SDS  digital  computer  in  accordance  with  operating  instruc- 
tions . 

2.  Light  off  ADAGE  graphics  computer  in  accordance  with  operating 
instructions  and  load  library  program  "Gated". 

3.  Load  the  program  in  the  XDS-9300  computer.   If  an  overlayed  version 
of  the  program  is  used,  the  entire  program  may  he  loaded.   If  an  overlayed 
version  of  the  program  is  not  used,  computer  memory  limitations  allow  loading 
only  two  antennas  at  a  time.   The  input  resistance  branches,  gain  branches, 
and  required  subroutines  for  the  antennas  desired  should  be  loaded  along 
with  the  main  program.   A  missing  lables  warning  will  result  but  the  program 
may  be  operated  if  only  antennas  loaded  are  called. 

4.  When  the  input  light  on  the  teletype  is  lighted  type  IDEV  =  1*  if 
ADAGE  1  is  to  be  used  or  IDEV  =  2*   if  ADAGE  2  is  to  be  used.   Pushing  the 
carriage  return  will  cause  the  data  input  format  to  be  displayed  at  the  graphics 
terminal. 

5.  Enter  parameters  and  option  commands  using  "Gated"  text  editing 
techniques.   Inputs  should  be  as  follows: 

a.  Under  ANTN  enter  one  of  the  following  to  specify  antenna  type: 

0001  Tilted  Dipole 

0002  Vertical  Whip 

0003  Vertical  Whip  with  Ground  Screen 

0004  Inverted  L 

0005  Sloping  Longwire 

0006  Sloping  Vee 

0007  Horizontal  Rhombic 

0008  Vertical  Half  Rhombic 

b.  Under  LENG  enter  length  in  format  F4.1 

c.  Under  HGHT  enter  height  in  format  F4.1 

d.  Under  PHIP  enter  4>  in  format  14 

e.  Under  THEP  enter  0   in  format  F4.Q 

f.  Under  FREQ  enter  f  in  format  F4.0 

g.  Under  EPSL  enter  e   in  format  F4.1 
&  r 
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h.   Under  SGMA  enter  a  in  format  F4 . 2 

i.   Under  PHI  enter  the  observation  azimuth  angle  for  the  vertical 
pattern  using  format  14. 

j .   Under  THET  enter  the  observation  zeniuth  angle  for  the  horizon- 
tal pattern  using  format  14. 

k.   Under  PARM  enter  0000.   If  reinitialization  is  desired  to 
erase  a  manually  entered  pattern,  enter  Q001.   If  Log  Gain  patterns  are 
desired,  enter  0002. 

1.   Under  ISTH  and  ISTV  enter  0000.  If  saving  the  pattern  that 

will  be  computed  in  the  current  compute  cycle  is  desired,  enter  0001.   If 

it  is  desired  to  keep  the  pattern  in  the  save  array,  these  option  commands 

must  be  set  to  0000  in  the  succeeding  compute  cycle. 

m.   Under  IRCL  enter  0000.   If  displaying  saved  patterns  is  desired, 
enter  0001. 

n.   Under  HGTT  enter  0000.   Entering  01.0  will  multiply  the  value 
of  sigma  by  .1.   Entering  02.0  will  multiply  the  value  of  sigma  by  .01. 

o.   Under  ALPH  enter  a  in  format  14. 

There  are  two  unused  data  blocks  which  no  operation  edits  must  be  made 
to  finish  the  data  input  processor. 

6.  Axes  and  a  blank  graphics  data  block  will  now  be  displayed  on  the 
terminal  screen.   A  pattern  desired  for  comparison  purposes  may  be  entered 
in  this  block  using  manual  graphics  editing  techniques.   To  terminate  this 
processor  operation,  push  the  end  edit  button  on  the  function  switch  panel. 
This  processor  may  be  terminated  without  entry  if  desired. 

7.  The  antenna  patterns  selected  will  be  computed  and  the  horizontal 
pattern  displayed  on  the  upper  axis.   Pushing  the  end  edit  button  will  cause 
the  vertical  pattern  to  be  displayed  on  the  lower  axis. 

8.  If  the  display  saved  patterns  option  has  been  selected,  pushing 
the  end  edit  button  two  additional  times  will  cause  the  vertical  and 
horizontal  patterns  to  be  superimposed  on  the  current  vertical  and  horizontal 
patterns.   The  program  will,  then,  return  to  the  enter  parameters  and  option 
commands  processor.   If  recall  has  not  been  selected,  the  program  will  return 
to  the   enter  parameters  processor  from  terminating  the  vertical  pattern 
display  processor  termination  (end  edit) . 

The  compute  cycle  is  now  repeated.   Ending  the  program  must  be  done  in 
accordance  with  laboratory  operating  instructions.   Figures  4.1  thru  4.20 
are  the  entries  for  the  examples  of  section  4. 

p.   Under  ISEA  enter  sea  state  in  14,  if  a  dynamic  display  is  desired 
for  dipole,  whip  or  longwire  antennas.   If  dynamic  display  is  not  desired, 
enter  0000. 

q.   Under  ICRS  enter  relative  direction  of  seas  if  dynamic  display 
is  desired. 
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